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We propose a variant of the CCS process algebra with new features aiming at allowing multiscale 
modelling of biological systems. In the usual semantics of process algebras for modelling biological 
systems actions are instantaneous. When different scale levels of biological systems are considered 
in a single model, one should take into account that actions at a level may take much more time 
than actions at a lower level. Moreover, it might happen that while a component is involved in one 
long lasting high level action, it is involved also in several faster lower level actions. Hence, we 
propose a process algebra with operations and with a semantics aimed at dealing with these aspects 
of multiscale modelling. We study behavioural equivalences for such an algebra and give some 
examples. 



1 Introduction 

Formal modelling notations of computer science are nowadays often applied to the description of bi- 
ological systems. Such notations can be used to unambiguously describe the structure and the events 
governing the dynamics of the systems of interest, thus allowing development of analysis tools such as 
simulators and of formal analysis techniques based, for instance, on model checking or on behavioural 
equivalences. 

As examples of formalisms that have been applied to the description of biological systems we 
mention Bio-PEPA (BlIISl, the (stochastic) TT-calculus ||29l [Ml |32l, Bioambients EH, the (C-calculus 
ifTTl . the Language for Biological Systems (LBS) (TJ), and the Calculus of Looping Sequences (CLS) 
ll22i r71[8ll. In these formalisms the dynamics of a biological system consists of a sequence of events (usu- 
ally biochemical reactions) described either as communications between processes of a process algebra, 
or as applications of some rewrite rules. In the stochastic extension of these formalisms, the dynamics 
of a system is described by taking also into account the different rates of occurrence of the events. Rates 
depend on some parameters associated with the events (e.g. kinetic constants of the corresponding bio- 
chemical reactions) and on the abundance (or concentration) of the entities (or reactants) that can cause 
such events. The rates are then used, as in Gillespie's algorithm lITSll . to describe both the exponentially 
distributed time elapsing between two subsequent events, and the probability of an event to occur. 
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This way of describing the dynamics of biological systems with sequences of events, however, as- 
sumes that the occurrence of one of such events can be described as an instantaneous change in the 
system state. In fact, even in the stochastic approach, the only notion of time that is considered is given 
by the frequency of the events, rather than by their duration. The duration of an event is usually ignored 
if it is negligible with respect to the time interval between two events, or hidden in such a time interval 
by chosing a rate for the event that is small enough to take into account both of its frequency and of its 
duration. 

Phenomena of interest in the study of biological systems often include processes at different levels 
of abstraction. A typical example is cell signalling, that involves gene regulation and protein interaction 
processes at the intra-cellular level, a signal diffusion process at the inter-cellular level, and some macro- 
scopic change at the tissue level. The processes at the different levels influence each other. However, they 
involve system components of very different sizes and are characterized by events having very different 
durations. This motivates the multiscale approach to modelling, whose application to biological systems 
seems to be promising Bll2ll9l [T3ll33l[34l . 

The description of the dynamics of biological systems by means of sequences of instantaneous events 
is not suitable for multiscale models. When different scale levels of biological systems are considered in 
a single model, one should take into account that events at a level may take much more time than events 
at a lower level. Moreover, it might happen that while a component (e.g. a cell) is involved in a long 
lasting high level event (e.g. mitosis), it is involved also in several faster lower level events (e.g. protein 
sinthesis) that neither consume such a component nor interrupt the higher level event. Consequently, the 
distinction between the time scales of the events, the possibility of having the same component involved 
in several events at different time scales, and the fact that the completion of some events may (or may 
not) interfere with other events in which the same component is involved require new notions of system 
dynamics to be considered. 

In this paper, we propose a process algebra with operations and with a semantics aimed at dealing 
with these aspects of multiscale modelling. Actually, we aim at undertaking a foundational study. Hence, 
we consider a minimal process algebra for the description of biological systems (a fragment of CCS [23l 
and of the Chemical Ground Form [14]) and we make the minimal changes we think to be necessary to 
describe the new aspects of interest. 

As regards the syntax of the new process algebra (called Process Algebra with Preemptive and Con- 
servative actions), we propose a new action prefixing operator that allows an action to be executed in a 
conservative (or non-consuming) manner, namely without removing the process that performed it. As 
regards the semantics, we define it as a labelled transition system by following the ST semantics approach 
(see ll2Tl[T9ll20l [T2'l) in which actions are not instantaneous, but described by two separate starting and 
ending transitions. This will permit to have processes in which multiple actions are running in parallel 
and competing for their completion. Indeed, we change the usual interpretation of the summation op- 
erator (by making it slightly similar to a parallel composition) order to allow a process to be involved 
in several actions at the same time. The termination of an action in a summation may interrupt (in a 
preemptive way) the others that are concurrently executed in the same summation, depending on which 
action prefixing operator is used. 

The semantics of the process algebra is given in a compositional way, and this allows us to study 
behavioural equivalences. In particular, we define a notion of bisimulation for the process algebra and 
we prove a congruence result for it. Some examples are given of use of the process algebra and of the 
bisimulation relation. 

The paper is structured as follows. In Section |2] we introduce the syntax of the processes, process 
configurations and some auxiliary functions to define the semantics, given in Section |3] with an exam- 
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pie. In Section |4] a behavioural equivalence for the processes is presented. Finally, we end with some 
conclusions and future work in Section [5] 

2 A Process Algebra with Preemptive and Conservative actions 

In this section we present the syntax of a Process Algebra with Preemptive and Conservative actions, in 
the following shortly denoted as PAPC. 

2.1 The Syntax 

PAPC is a process algebra with dyadic communication in the style of CCS |[23l and of the Ti-calculus 
Il24ll25]| . Hence, communications involve exactly two processes at a time. As in the cases of CCS and 
of the ;r-calculus, the description of biological systems with PAPC is based on a processes-as-molecules 
view: the modeling of a molecular species S occurring n times in a system includes n copies of a process 
Ps modeling a single molecule. 

Differently with respect to the approach of classical process algebras, where actions are instanta- 
neous, in PAPC we assume that actions can consume time and that the instants of start and completion of 
an action can be detached. An action that has already started but not yet completed is said to be running. 

Let us assume an infinite set of actions Act ranged over by a, j8, . . . and a function (-) : Act Act 
such that a = a. We also denote with Ac?^ the set of actions enriched with the special internal action T, 
Act-c = Act U {t} . The abstract syntax of PAPC is a follows. 

Definition Processes of PAPC are defined by the following grammar: 

P := I a.P I a:P \ P + P \ P\P \ A 

where a G Act. We denote the set of all processes by 

As usual, denotes the classical idle process that can perform no action. PAPC processes can perform 
actions mAct in two different ways (represented by two different action prefixing operators). Both a.P 
and a : P can perform the action a. However, after performing such an action, the first process simply 
behaves as P, while the second process continues its execution as P | a : P, namely it produces a copy 
of P, but it is also availabe again to perform a. This difference in the behaviour allows a.P to model a 
molecule that can be involved in a reaction that transforms it into another molecule, and a : P to model 
a molecule that can be involved in a reaction that does not consume nor transform it. Another difference 
between a.P and a : P is related with the fact that actions are not instantaneous and that a process can 
start several actions concurrently. When an action used as in a.P completes, it interrupts all other running 
actions for the same process. This agrees with the intuition that a.P represents a transformation of the 
described molecule into a different one. On the other hand, when an action used as in a : P completes, it 
does not interrupt the other running actions of the same process. Again, this agrees with the intuition that 
the described molecule is neither consumed nor transformed into something different. We say that a is 
preemptive with respect to a.P and that it is conservative with respect to a : P. Notice that in PAPC the 
composition of these two approaches is possible, since an action a may be conservative for a process, 
and the complementary action a may be preemptive. 

A process P + 2 is able to start actions of both P or Q, but notice that, because of the difference 
between conservative and preemptive actions and the fact that actions are not instantaneous, this summa- 
tion operator does not correspond to the choice operator of classical process algebras. In fact, starting an 
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action in P (resp. Q) does not imply that process Q (resp. P) is discarded. As a consequence, P + Q may 
start several actions, which are said to be in competition, meaning that they run concurrently until one of 
them completes. When an action a in P (resp. Q) completes, then all actions in Q (resp. P) competing 
with it have to be interrupted only if a is preemptive. If this is the case, then Q (resp. P) is discarded. 
Differently, if a conservative action completes, then the summation is not discharged and, consequently, 
all the running actions in the summation are still running. 

A process P \ Qis the parallel composition of processes P and Q. Handshaking is possible between 
any action a in f and its complementary action a in Q. The handshaking is to be performed to assign to 
an instance of communication a unique identifier which may be used to interrupt the communication. An 
interesting case here is when an action a in P and its complementary action a in Q have been coupled 
with an handshaking, have started and not yet terminated, and some preemptive action in P competing 
with a completes, so that a must be interrupted. In fact, in this case also the complementary action a in 
Q must be interrupted. 

Finally, constants A are used to specify recursive systems. In general, systems are specified as a set of 

def 

constant defining equations of the form A = P. As usual, we assume that all processes in these equations 
are closed and guarded. 

Some further considerations are worth in order to introduce the differences between PAPC and the 
classical process algebras. The capability of having competing actions is at the basis for the choices we 

def 

made in the definition of PAPC. A process P = Pi + ... +P„ can start multiple actions in parallel, but can 
be involved in each action at most once at a time. Notice that, in classical process algebras, this is not 
possible since an action, when starts, determines the future process to transform P in. Hence, the choice 
is resolved at the time of the starting of an action. In this sense, the summation operator of PAPC is not 
a classical choice for the reason that the competing actions compete for their completion, and, then, the 
semantics of the completion, and hence the semantics of the PAPC summation, will depend on the type 
of the action to be completed, namely whether it is conservative or preemptive. 

Consequently, at any time of a computation, P could be in a configuration in which some of its actions 
are currently running. More precisely, the competition of the running actions is due to the fact that they 
are waiting to complete. Practically, the time for completion may be modeled by general distributions 
as in [[12], or by delays as in [5, 6J. However, in this first definition of the algebra we do not consider 
quantitative timing and stochasticity. 

In order to define the semantics of PAPC we need to model a process with possibly running actions. 
We do this by introducing a notion of process configuration. 

Definition Process configurations of PAPC are defined by the following grammar: 

Cp := [a]' .P I [a]':P | Cp + Cp \ Cp\Cp \ P 

where a G Act and / G N. We denote the set of all possible process configurations as 

Any process P G ^ is also in a valid configuration, hence C'^. However, a process configuration 
may contain actions denoted by a different prefix. In particular, the configuration [a]' .P is the config- 
uration reached by a.P after a has started, and [a]' : P is the configuration reached by a : P after a 
has started. For both the action prefixes, the new argument Z G N is a natural number that identifies the 
running action. Notice that these identifiers, which have to be unique, are computed by the handshaking 
performed before the start of an action and, once a preemptive action is completed, they may be used 
to interrupt all other competing actions. Then, if one of these competing actions is a and it has started 
an handshaking with another action a in a process running in parallel, then also this action a will be 
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interrupted. This can be obtained by assigning the same identifier to these two actions when the hand- 
shaking begins. By the definition of the semantics it will be clear how both the partners will share the 
same identifier for the actions. 

Let us define by structural recursion an auxiUary function W : H- N as follows: 

Id{[a]' .P)=Id{[a]' :P) = {1} 
ld{Cp + C'p)=Id{Cp \ Cp)=Id{Cp)Uld{Cp) 
Id{P)=<d. 

The value ld{Cp) denotes the set of the identifiers of the actions in the configuration Cp that are running. 
For instance, given a configuration Cp == [a]' .P + j3 : 2| [7]' .T , the identifiers collected by function Id 
are given by Id{Cp) = {1,1'}. 

Finally, we define a function Action : N x i-)- p[Act) such that Action{l,Cp) collects the set of 
actions currently running in Cp and with assigned identifier /, if any. The function Action is defined as 
follows: 

Action{l, [a]' .P) =Action{l, [a]' : P) = {«} 

Action{l, [a]'' .P) = Action{l, [a]'' ■.P)=% if / 7^ /' 
Actional, Cp + C'p) =Action{l,Cp \ Cp) = Action{l,Cp)U Actional, Cp) 

Action{l,P) =0. 

The definition of this function is similar to the definition of Id. For instance, given the configuration 
Cp == [a]' .P + j8 : Q\ [7]' .T, the actions collected by function Action are given by Action{l,Cp) = {a}, 
Action{l' ,Cp) = {7} and, for all /" / / and /" / Action{l" ,Cp) = 0. 

In the next sections we define the semantics of PAPC by using the notions of process, process con- 
figurations and these auxiliary functions. 



3 A Structural Operational Semantics for PAPC 

In this section we define a Structural Operational Semantics (SOS) ||28l for PAPC. The aim of the SOS is 
to equip PAPC with a Labeled Transition System (LTS), namely a set of transitions of the form P \y P' 
representing a move from P € to G with the label £ carrying some information about the move 
and the index r used to group transitions describing a particular aspect of the behavior of the processes. 
The LTS is defined by a set of SOS transition rules of the form 3^^^^- Intuitively, each of these rules 
explains how a move of a process is obtained from moves of its subprocesses. All our rules are in Figures 
[T]-[6] We assume the standard way for assigning an LTS with such a set of transition rules (see, e.g., 111). 

The main features of the SOS we want are the following. Firstly, it must have a mechanism to 
interrupt competing actions and this mechanism is activated by the completion of a preemptive action. 
Secondly, the style of the semantics must be ST-like, as this permits to easily observe detached events as 
the start and the completion of an action. 

In order to get this features, we define a relation for modeling the start of an action and the coupling of 
processes; this will be named as the handshaking relation. Furthermore, we define a completion relation 
for modeling the finishing of both preemptive and conservative actions. These two relations will make 
use of an interruption relation to model the interruption of currently running actions, as required by the 
notion of preemptive actions. 
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(H,) a.P^H[a]'.P {H2) a:P^H[a]' -.P 



P^hP' HldjQ) 
P + Q^hP' + Q 



(H 



P^hP' leId{Q) V = min{n-Id{P + Q)} 
P + Q^hP'[1'/1] + Q 



(H5) 



P^hP' HId{Q) a^AcU 

p\Q^hP'\Q 



P^hP' leId{Q) I' = min{n-Id{P\Q)} aeAct^ 



P\Q^hP'[1'/1]\Q 
P^hP' Q^hQ' I" = min{n-Id{P\Q)} 

p\Q^HPT/i]\QT/n 



Figure 1: The handshaking relation — x 0+ x 
The handshaking relation 

This relation is used to model the starting of an action and the coupling of the processes starting comple- 
mentary actions. The handshaking relation is — >//C 'rf x 0+ x 'rf, where 0+ contains labels 0+ of the 
form 

0+ = (/,«+) 

where Z € N represents the identifier assigned to the started action a G Act^, and the use of the superscript 
"+" comes from the definition of the semantics in the ST style, in order to denote the start of an action. 
The SOS rules in Figure [T] are at the basis of the definition of — We implicitly assume the rules 
symmetric to (//g), (7/4), (^5), (^e)- 

Rules (Hi) and {H2) model the starting of an action a. At any time a process with prefix a can 
start action a moving to a configuration in which it cannot perform the same action anymore, i.e. the 
configuration [a]Vp or, analogously, the configuration [a]' : P. Such a configuration, together with 
the one describing the process performing the complementary action, has to be uniquely identified by 
a natural number representing the identifier of the just started action. At this step, the process simply 
chooses 1 as unique identifier. All our choices for assigning identifiers to actions are inspired to those 
of [11], which ensure that the portion of LTS rooted in a given process is finite. The rules for binary 
operators + and | will solve conflicts of colliding identifiers, if any. Notice that both preemptive actions 
and conservative actions start in the same way. 

Rules (7/3 ) and {H4) combine the start of an action with operator +. In rule (7/3 ) the identifier I of 
the action a started by process P has no conflicts with the identifiers of the competing actions running in 
process Q. Differently, in the case of rule {H4) a conflict does exist, which implies that a fresh identifier 
/' replaces /. Again, the policy by which we choose the new fresh identifier, along the line of [1 1 1, is such 
that the resulting LTS is finite. More precisely, the use of the set N — Id{P + Q) is such that we consider, 
in the process of renaming an identifier in P', the only set of identifiers not used in a process P + Q and. 
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(/i) [a]'.P^ia.P ih) [a\':P^ia:P 
[h) [atP%[a\'.P ih) [a]' ■.p\i[at.P 



{I5) a.P^ia.P (k) a:P^ia:P 
P + Q^iP' + Q' 



Figure 2: The interruption relation — >/C x p{N) x ^. 



from that, the choice of extracting the minimum value is such that the LTS is finite. We use this strategy 
in all the rules where we have to resolve some conflicts. 

Rules {H=,), {H(,) and {Hj) combine the start of an action with the operator |. Rules {H5) and (//g) 
model an autonomous move by one of the two processes, and deal with identifiers as {Hj) and {H4), 
respectively. As in classical process algebras, we do not force P and Q to handshake, since P could 
handshake with a further process composed in parallel with P \ Q. 

Notice that here we may have a conflict even if in Q the action associated with the colliding identifier 
is the complementary action a. Rule {Hj) models the handshaking by assigning to this particular instance 
of synchronization a new fresh identifier /" chosen with the same policy used to resolve conflicts in the 
previous rules. The renaming of both old identifiers with the newly generated is due to the fact that, in 
general, the two processes will have two different candidate identifiers, i.e. / and The system in this 
case exhibits the internal action T+. By applying this rule, the two processes terminated this handshaking 
phase. 



The interruption relation 

This relation, differently from those found in classical process algebras, is used to model the interruption 
of a set of actions currently running in a process. The interruption is caused by the completion of 
competing preemptive actions. The interruption relation is — >/C 'if x p(N) x 'rf, where a label M € p(N) 
contains the identifiers of the actions that have been interrupted. The rules presented in Figure [2] are at 
the basis of the definition of — >/. 

At any time, a process either in configuration [af .P or [a]' . P may interrupt the action it is currently 
performing. In these cases, treated with rules (/i) and (72), it moves to a configuration in which the 
interrupted action a may start again, namely to configuration a.P or a . P, respectively. In both the rules, 
the identifier / of the interrupted action is exhibited as a label of this transition. Again, this information 
will be used to interrupt also the partner of this action, as we are assuming that there is a partner in the 
system which, after terminating the handshaking phase, has been coupled with the same label /. 

In some cases not all the actions have to be interrupted, so the processes in configuration [a]' .P or 
[a]' : P must be able also to non-deterministically decide whether to interrupt or not. This second case is 
described by rules (It,) and (74), which may seem controversial at first glance. In particular, it may not be 
clear why a process may independently decide whether to interrupt or not some of the currently running 
actions. The need of this autonomy for the process can be clarified by an example. Let us assume a 
process configuration {P + L) \ {Q + \ S \ R, where both P and Q successfully complete a preemptive 
action. The actions to be interrupted are those currently running in both £ and namely those with 
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(CO [a]'.P^cpP (Q) — 



p+Q !^^^:^^cpP' 



.^.P^^CpP' Q^iQ' MD{Id{Q)nL) V = {LUM)\{Lnld{Q)) a€Act, 

^^^> ] — \ 

P\Q^~^cpP'\Q' 

. . P^cpP' Q^^cpQ' N = LnM 

PIQ^^^^^^^CPP'IQ' 

Figure 3: The completion relation for preemptive actions ->cpQ ^ x ®cp x ^■ 

identifiers denoted by Id{L) Uld{L'). Let us assume that some of the actions that have to be interrupted 
in E and Z' were coupled with some actions in S. In this case, also these actions in S should be interrupted 
as well. Moreover, S may be involved in other actions currently running and coupled with actions in R. 
Indeed, these actions must not be interrupted. This means that from S the correct derivation with the 
interruption relation, in general, will not exhibit as label Id{S), indeed it will exhibit a strict subset of 
Id{S). This implies that S must be able to autonomously decide which actions to interrupt, and this can 
be done by properly combining derivations of the interruption relation. The composition of the relations 
of the whole semantics will provide the correctness, namely the fact that all and only those to interrupt 
are actually interrupted. 

Also, a process which is not performing any action, namely a process in a configuration a.P or a : P, 
does not interrupt any action, as stated by rules (/j) and (/g). 

Finally, rule (/y) simply collects the labels of the interrupted actions in a summation. Notice that this 
relation is not defined for process configurations of the form Pc \ P'q as the use of this relation is limited 
to the level of the summation. 



The completion relation for preemptive actions 

This relation is used to model the completion of a preemptive action. We will define completion relations 
also for conservative actions as well as the combination of both preemptive and conservative actions. 

The completion relation for preemptive actions is — >cp^ x ®cp ^ "^^^^ ®cp containing labels 
of the form 

e ■= {l,a-,N) 

where Z G N represents the identifier that was assigned to the completed action a € Acf^ when it was 
started, N € p(N) is the set of the identifiers of the competing actions that are interrupted by the termi- 
nation of a, and the use of the superscript "-" comes from the definition of the semantics in the ST style. 
The rules presented in Figure [3] are at the basis of the definition of — >cp- We implicitly assume rules 
symmetric to (C2) and (C3). 

Rule (Ci) describes the completion of a preemptive action. When it completes, as the action is 
preemptive, the process is substituted by its continuation P. In the label, the identifier / is needed to 
couple this process with the one performing the corresponding complementary action a, which will have 
the same identifier / because of the handshaking, and states that no action is interrupted. 
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(C5) [a]':P^-^^^cc«:P (C, 



6) 



I LUM P" 

P + Q ' ' ) ccP' + Q' 



... P^^^ccP' ^^. P^^ccP' q'-^^ccQ' 

^^^> !a-LP" 77~0 

P I Q ^'^^cc P'\Q P\Q -^cp P'\Q'\ P" I Q" 

Figure 4: The completion relation for conservative actions — >cc^ x ®cc ^ ^• 

Rule (C2) states that the completion of a preemptive action in P affects a summation P + 2 so that 
all actions running in Q should be interrupted. This is obtained by adding to the set of labels of actions 
interrupted L, the set of actions currently running in the process Q which disappears by the completion 
of the action in P, hence the exhibited set of labels becomes LUld{Q). 

Rule (C3) states that the completion of a preemptive action in P affects a parallel composition P \ Q 
so that all actions running in Q that are coupled with actions interrupted in P, must be interrupted as well. 

Rule (C4) models the case in which both P and Q complete preemptive actions that were coupled. As 
in classical process algebras, the whole system P | Q exhibits an internal action T. Some of the actions 
required to be interrupted outside P, may be also required to be interrupted by Q. Such a set is denoted 
by A'^ and can be removed from the set of actions that can be interrupted outside P \ Q. The remaining 
set of actions, which have to be still interrupted by further composition with the parallel operator outside 
P\Q, is the set of those belonging to P and not to Q, and viceversa. 



The completion relation for conservative actions 

This relation is used to model the completion of a conservative action. This relation is — >cc^ x ©cc ^ 
where 0^(-. contains labels of the form 

e := {l,a-,N,P) 

where / € N represents the identifier assigned to the completed action a € Acf^, A'^ € is the set of 
identifiers of the interrupted actions, and P ^ is the continuation of the action which terminated and 
that, syntactically, must be propagated at the level of a parallel composition. At first sight it could sound 
strange that we need the component N. The idea is that we have to take care that a process Q terminating 
a conservative action a could be composed in parallel with another process Q' terminating the action 
a coupled with a. Now, if a is preemptive, there may be some running actions j8 in Q' that should be 
interrupted, which implies that if there is an action j8 in 2 coupled with j8 , also jS must be interrupted. 
For this reason, such a p must appear in N. Of course, if also a is conservative, then in the transition by 
Q used to infer the transition of Q\ Q', N will be empty. 

The rules presented in Figure |4] are at the basis of the definition of relation — >cp- We implicitly 
assume rules symmetric to (Cg) and (C7). 

Rule (C5) deals with termination of a conservative action a in configuration [a]' : P. The process 
becomes able to perform the terminated action again, namely it rolls back to configuration a : P. Also, 
as expected by a conservative action, it produces the continuation P, which, because of the inductive 
approach of the SOS semantics, cannot appear at the same syntactic level of the configuration a : P. 
Specifically, the continuation P will have to appear at the level of a parallel composition. To forward P 
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.^. P^^^ccP' Q^CPQ' L<ZM 
p\q'-^cpP'\Q'\p" 

Figure 5: The completion relation for hybrid actions obtained by means of the other completion relations. 

AArP' 

Figure 6: The standard rule for recursion. 

at the coiTcct syntactic level, P is exhibited as a label of the transition. The empty set used in the label 
denotes that no action is interrupted by completion of a. 

Rule (Ce) clearly justifies the terminology "conservative". When a conservative action a completes 
in a process P being part of a configuration P + Q, it is required neither that actions in Q are interrupted, 
nor that Q is canceled. This is clearly different from what happens when a preemptive action is completed 
(rule (C2)). More precisely, the set of actions interrupted in Q, namely M, is a subset of Id{Q) since here 
not all the actions in Q have to be interrupted. Furthermore, the continuation of the action, namely P" in 
the rule, is exhibited as a transition label, since also in this case we are not yet at the syntactic level of a 
parallel composition. 

Rule (C7) describes the case in which P completes a conservative action a and the coupled action a 
is not in Q, since it runs in some other process running in parallel with P \ Q. 

Rule (Cg) deals with the completion of two coupled conservative actions. As expected, the system 
exhibits an internal action T, the label shows that no action has to be interrupted, and both the continua- 
tions appearing in the labels of the transitions of both processes P and Q, namely P" and Q", are put in 
parallel with the continuations P' and Q'. Notice that this last derivation is a derivation for — >cp rather 
than — >cc- The reason for this is that P" and Q" are already at the correct syntactic level and do not 
require to be lifted anymore. 

Completion of both conservative and preemptive actions 

We have to deal with the completion of two coupled actions a and a such that one of them is conservative 
and the other preemptive. To this purpose, we add the rule shown in Figure |5] and we implicitly assume 
a symmetric rule. 

Notice that L QM expresses that the actions running in the process Q performing the conservative 
actions and that have to be interrupted are those that were coupled with actions running in the process P 
performing the preemptive action. In fact, such a coupling is the only reason we have to interrupt actions 
running in Q. 

Recursive definitions 

As far as the naming of processes is concerned, we define the standard rule for recursion, showed in 
Figure |6l for all the relations we defined. 
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A toy example 

In this section we discuss the modeUng of a multiscale system where we consider two populations. At 
a higher level of abstraction we consider a cell C, and at a lower level a generic protein P. A cell C can 
be involved in a process leading to its duplication. Also, it can be involved in some low-level reactions 
(i.e. DNA transcription inside its nucleus) leading to the creation, in the environment outside C, of a 
protein of species P. Of course, we consider this model at a level of detail such that we do not need to 
take into account any other possible population of either cells or proteins which could be involved in the 
dynamics. 

In the context of chemically reacting systems such a system may be described by two populations C 
and P, and by two reactions Ri and R2 such that 

Ri:C^C + C R2:C-^C + P. 

Reactions R\ and R2 model the non linear growth of cell C and the production of a protein P by a cell C, 
respectively. Notice that, at this level of detail, the production of protein P depends on the cell C where 
all the details of the biological process leading to the creation of the protein are abstracted away. The 
initial state of the system can be defined to have a precise initial number of cells C and proteins P in the 
environment. 

We model now such a system in PAPC, and we show how the semantics models the behavior of the 
populations. As PAPC is based on the paradigm processes-as-molecules, we start by assuming two types 
of processes for each species which, for clarity, are named C and P. As we want to model two reactions, 
we assume the following set of actions {a, a, 7,7} where a (resp. a) and 7 (resp. 7) model reaction 
Ri and R2, respectively. Also, as in PAPC the communication is dyadic and the reactions use a single 
reactant, we define two auxiliary process X" and X"^, used to model the communications on a and 7, 
respectively. 

Reaction Ri creates two new different cells able to start again, if possible, the duplication process. 
In the context of PAPC we model Ri by using preemptive actions for both a and a since the duplication 
of a cell interrupts all the low-level protein-transcription event inside the duplicated cell. Differently, the 
actions modeling reaction R2 are conservative since the protein-transcription event does not interrupt the 
duplication process started by a cell. 

The PAPC processes are defined as 

C=a.{C\C) + Y:P X« = a.(X«|X«) X^=y:0. 

Notice that we do not give a definition of the process P since it appears only as a product of the events 
we want to model, and hence we are not interested in the interactions it may have in the system. 

Process C as expected can perform two actions. Action a produces the two new copies of C. Such 
an action is performed by synchronizing with the auxiliary process X" which will produce two copies of 
itself to permit the duplication of the new cell. Indeed, the number of copies of the auxiUary process X" 
must grow with the same law of growth for the cells C. Both C and behave as preemptive in a and 
a, as expected. Process C can also perform, by synchronizing with X'^, the action 7. The result of such 
action is to produce a new protein B without interrupting its running duplication, if any. The number of 
processes in the system bounds the number of cells which can simultaneously produce a protein of 
type B. This is obtained by producing, in X'^, the nil process. If this were a too strong constraint, it would 
have been possible to add a process in the continuation of the action a in C, as done for X". For the 
sake of simpUcity in this model we consider this constraint to be reasonable. 
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We discuss now some features of the semantics of PAPC for a simple system 5 described by the 
process 

In S both the reactions may fire. We assume reaction Ri to fire first. To this extent, the semantics permits 
to observe handshaking derivations as 

C^H[a\\{C\C) + Y:P Z« | X«) . 

Then the whole process 5 performs a derivation as 

5-^J^[a]^(C|C) + 7:P| [a]V(Z" |Z^ = 5' 

where the new process 5' is such that the action a is now running in C and in X", with identifier 1. In 5' 
action a can not start, but just complete, we assume reaction R2 to fire. The semantics permits to observe 
the handshaking derivations 

The composition of these derivations resolves the conflicts of the colliding identifiers such that the deriva- 
tion for C will be 

[a]\{C\C) + y:P^H [a\\{C \C) + [yf -.P 
and the whole system performs the following derivation 

S' [a]' .(C I C) + [yf : P \ [a]' | Z«) | : = 5" 

where in S" all the possible actions are running. We consider now two different cases: (a) R\ completes 
before R2 and {b) viceversa. 

(a) Reaction /?i completes before R2. in this case action a (resp. a) completes before action 7(resp. 
7), interrupting it. The semantics permits to derive transitions as 

[af .(C I C) ^^cp C\C [a]' .(^" I X") ii^cp^" I ■ 
Actions in C and in have to be interrupted, hence we derive 

[«]' -iC I C) + [rf : P -^^^-^cp C I C [7]' : ^/X^ 

where {2} denotes the actions to be interrupted. Consequently, the whole process 5" will perform 
the transition 

S" ^^cpC\C\X''\X''\X^ 

where in the resulting process no actions are running, as expected, and there are two cells and two 
auxiliary processes X". 
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{b) Reaction R2 completes before Ri : in this case action 7 (resp. 7) completes before action a (resp. 
a). The semantics permits to derive transitions as 

[7f:0™>cc7:0 [y]- : P '-^cc 7 : P . 

No actions have to be interrupted in any process, hence we derive 

[a]' .(C I C) % [a]' .(C I C) [a]' .{C\C) + [yf:P ^^^cc [a]\{C \C) + y: P . 

The whole process 5"' will then perform the transition 

S" ^^cp [a]^ .{C\C) + y:P \ [a]^(X« | X«) | | P [ 

where, as expected, in the resulting process only one action is still running (cell division) and a 
single protein P has been produced. 



4 Bisimulation equivalence for PAPC 

Bisimulation equivalence is a central notion in concurrency theory. For processes with a higher order 
behaviour, namely processes whose behaviour is described by portions of transition systems in which 
processes appear in the labels, the notion of bisimulation is usually replaced by a higher order notion of 
bisimulation |[3l|T0l|35l, which can be rephrased in our setting as follows: 

Definition A symmetric relation =^ C x is a bisimulation iff whenever {P, Q) € then it holds 
that: 

• HP 4, P' for any e 'T, r G {H,I,CP} and label £, then Q 4, Q' for some Q' e'^ such that 
{P',Q')e^. 

• ifP '^''"') cc P' for any P' € ^, / € N, a e Act, L C N and P" G ^, then Q '^'^") cc O! for 
some 2', 2" G such that (P', 2') G ^ and (P", 2") G ^. 

The union of all bisimulations is, in turn, a bisimulation, which is denoted with k, and is called as "the 
bisimulation". 

For an algebric treatment of bisimulation equivalence and to reason in a compositional way, a bisim- 
ulation is required to be a congruence. By taking the standard notion of context C[ ], a bisimulation M is 
a congruence with respect to all operations of the process algebra if and only if, given any pair (P, 2) G ^ 
and any context C[ ], it holds that (C[P],C[P]) G M. 

Theorem Bisimulation is a congruence w.r.t. all PAPC operations. 

Proof. If we consider PAPC without recursion, then the proof comes for free. In fact, in f26l it is shown 
that higher order bisimulation is a congruence with respect to all process algebra operations whose se- 
mantics is defined through transition rules respecting some syntactical constraints, and it can be checked 
that the transition rules we use respect such constraints. The extension of the proof to the case of recur- 
sion is standard. □ 
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The bisimulation relation for PAPC is a very fine beliavioral equivalence. Tliis can be seen as a dis- 
advantage, since with behavioral equivalences it is often desirable to be able to equate as many processes 
as possible. On the other hand, the fact that bisimulation turns out to be fine may have another meaning, 
namely that all the ingredients used in the process algebra play an important role. This does not hap- 
pen, for instance, for the parallel composition in some variants of CCS where it can be reduced into an 
equivalent summation of processes. 

We go through this last consideration via an example. Let us consider the following two PAPC 
processes 

C^^a.{C\C) C'^^a-.C. 

The behaviour of the two processes is similar: both of them perform an action a and then continue 
as with copies of the initial process (in the case of C' one of such two copies is not indicated in the 
continuation of a since it is obtained from the semantics of conservative actions). Even if the behaviour 
of the two processes seems to be the same, it is immediate to see that the two processes are not bisimilar. 

In fact, they repeatedly perform the same handshaking transition - — >h, but followed by two different 

completion transitions, namely ->cp and )-cc- This permits to state that, in general, C 56 C . 

This is exactly what we expect from our bisimulation relation since, when the two processes are put 
in a summation context, their behavior would determine the behavior of the whole context. In fact, 
the completion of the action in a process C + £ would interrupt any action currently running in £ and, 
differently, for the case of C' + £ no actions in £ would be interrupted. This is due to the fact that the two 
processes perform the same action but with different prefix operators. 

In order to see whether this is important, is enough to consider the toy example given in the previous 

def 

section where actions are modeled by a process C = a.{C \C) + Y '■ P- In that example the creation of 
two new cells should be an event which interrupts, when completes, the production of protein P by the 

def 

cell. Of course, if the process used in the model would have been C = a.C + y.P, then the completion 
of the duplication process for the cell would not have interrupted the production of P. 



5 Conclusions 

In this paper we considered the problem of modeling biological systems in which different scale levels 
are taken into account resulting in the fact that actions at a higher level may take more time than actions 
at a lower one. 

In order to model such systems we defined PAPC, a variant of the CCS process algebra in which a 
process may be simultaneously involved in one long lasting high level action and in several faster lower 
level actions. In order to model this, we added in the algebra two different prefix operators to model 
the role of a process in an action. A process can either act as conservative or preemptive in an action, 
resulting in two different behaviors for the process and the other actions which currently are started and 
not completed. 

We gave a compositional Structural Operations Semantics for PAPC by means of different relations, 
one for each of the possible events which may change the state of a process, namely the start and the 
completion of the actions. The semantics we gave is in ST style as this permits to observe the start and 
the completion of an action as two detached events. This style of the semantics permits also to observe 
processes in configurations in which multiple actions are started and not completed. 
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We also defined a notion of behavioral equivalence for PAPC processes based on the ideas of higher- 
order bisimulations for process calculi. We proved that our bisimulation is a congruence for all PAPC 
operators. 

In the paper, we also showed some simple example of PAPC processes such that their semantics 
permits to observe the key features of the algebra. We also discussed the notion of bisimulation we 
introduced by analyzing two simple PAPC processes. 

As a future work we will apply PAPC to the modeling of multiscale systems in order to prove the 
utihty of the formahsm. Also, we may consider to enrich PAPC with biologically inspired operators 
to easily model complexation, de-complexation or more complex biological structures as membranes 
or compartments as it has been previously done with other calculi. Moreover, we may consider the 
definition of more biologically inspired notions of equivalence for PAPC processes. 
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